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Abstract 

In the present article we investigate the influence of the contact region on the distribution of the 
chemical potential in integer quantum Hall samples, as well as the longitudinal and Hall resistance 
as a function of the magnetic field. First we use a standard quantum Hall sample geometry and 
analyse the influence of the length of the leads where current enters/leaves the sample and the ratio 
of the contact width to the width of these leads. Furthermore we investigate potential barriers 
in the current injecting leads and the measurement arms in order to simulate non-ideal contacts. 
Second we simulate nonlocal quantum Hall samples with applied gating voltage at the metallic 
contacts. For such samples it has been found experimentally that both the longitudinal and Hall 
resistance as a function of the magnetic field can change significantly. Using the nonequilibrium 
network model we are able to reproduce most qualitative features of the experiments. 
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I. INTRODUCTION 



Nowadays the steps in the transversal (Hall) resistance at the inverse of integer multiples 
of e 2 /h are well understood in terms of one-dimensional edge channels The number 
of such channels decreases with increasing magnetic field due to increasing energy spacing 
between the Landau levels. 

Despite this clear physical picture it was found experimentally that width and magnetic 
field values of the transition region between plateaus can change in case of nonideal contacts, 
e.g. when applying gating to the contacts. Since a few years experiments using scanning 
probe techniques (see Ahlswede et. al. [2|) were performed to visualize the distribution of the 
nonequilibrium potential near non-ideal contacts. On the other hand also the distribution 



in order to optimize the contact 
I) found a strong indication that 



of current in quantum Hall samples was investigated |3[ 
geometry. A significant number of experiments (see |4|-[( 
the contacts do not behave as ideal ones, meaning that not all edge channels seem to arrive 
at the reservoir of the metallic contacts of the sample. In addition imperfect equilibration 
among edge channels can also lead to deviations in the measured resistance values. Fur- 
thermore, deviations of the R xx peaks from the expected shapes were observed for non-ideal 
contacts by Dahlem et. al. 6], who proposed that the deviation might result from the crys- 
talline orientation of the edge of the contacts. The dependence of the contact resistance 
on the crystal orientation has also been subject to other experiments They embedded 
Au/Ge/Ni contacts in various Alo.Q7Gao.33 As /Ga As heterostructures and found that the 
contact resistance varied with length and orientation of the interface line of the contact. 

The experimental findings show nontrivial results for the potential distribution and resis- 
tances as a function of details of the sample geometry. In order to simulate such samples we 
use the nonequilibrium network model (NNM) (Ref. [8j and 9j) because the actual geometry 
of the sample can be taken into account. In addition the NNM has already proven successful 
to simulate even complicated sample geometries like anti-Hall bars 101 ] . 

Regarding theoretical investigations of nonideal contacts in simplified models, van Wees 
et. al. discuss influences of non-ideal contacts in terms of a nonzero reflection probability of 



electrons at the contact, which leads to waves traveling in the system This leads to 
potential resonances and a fine-structure of the conductance between plateaus. We did not 
consider such effects here and assumed zero reflection at all our contacts. Note in case of 
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slowly varying electrostatic potential from the system into the contact the reflection would 



12|. 



be suppressed 

The article is structured as follows: In Sec. [IT] we present the setup and underlying theory 
of the NNM. We investigate the influence of the geometric shape of the contact and the 
current-injecting leads as well as potential barriers within leads and measurement contacts 
in Sec. IIII Al of the results. Simulations of gated samples are then presented in section UlI B[ 
Finally we summarize and conclude in Sec. IIVI 

II. THEORY 

The exact Hamiltonian of an integer quantum Hall sample is not easily formulated be- 
cause the disorder potential (long ranged), which defines the underlying electrical potential 
besides the atomic (microscopic) environment and the confining potential (due to interfaces), 
depends on details of the exact configuration within the sample. This lead to a variety of 
more or less simplifying models. The most prominent of this type of models is the Chalker- 
Coddington model [13] , which is able to predict delocalized states but cannot describe the 
nonequilibrium steady state. However it seems to lead to the correct exponent of the corre- 



lation length [14H16| from which we conclude that the topology of the setup of the network 
is sensible even in the nonequilibrium regime. 

Due to the strong magnetic fields Landau levels develop with localization within the mag- 



netic length Ib = \fhfeB given by the cyclotron motion, e and B denote elementary charge 
and magnetic induction. Considering the very different length scales involved, electrons 



are confined to move semi-classically along equipotential lines of the disorder potential 17 ]. 



Regarding this approximation transport occurs via tunneling constrictions, that is, saddle 
points in the underlying potential landscape. The elastic tunneling transition probabilities 
across such saddles have been calculated quantum mechanically for equilibrium 18, Ji]]. 

In order to be able to discuss special wiring configurations we have to define the notion of 
longitudinal and transversal components with the help of a curvilinear coordinate system, 
determined by the local orientation of the contour of the chemical potential in the plateau 
regime. In this way the local longitudinal direction is defined within the plateau state 
via the local tangent to the equipotential line. Clearly no distinction between longitudinal 
and transversal resistance can be made in the transition region. If then two measurement 
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contacts can be connected within the plateau region along an equipotential line, we measure 
a longitudinal resistance. In case of intersections a transversal (Hall) component is measured. 
Note however, the NNM calculates within a fixed coordinate system labeled by (x,y). As 
required the resistances do not depend on the choice of these coordinates. 

The physical content of the NNM can be understood in terms of local equilibrium (see 
e.g. Ref. [20] for a formulation in continuous space variables) when reformulating it for the 
network of semiclassical wavefunctions. In this way we attribute unique thermodynamical 
quantities such as the chemical potential to each single wavefunction, given by a trajectory 
along the contour connecting two saddle points. 

Oh and Gerhardts 21] use a very related concept of local equilibrium, however in a contin- 



uum description, together with the Thomas- Fermi or Hartree approach to the self-consistent 
calculation of charges. They use a simple model geometry and assume translational invari- 
ance in longitudinal direction which results by continuity in a current independent of the 
transversal direction. They also include the velocity of particles when averaging over states 
to obtain the charge density. Current injection and the important influence of contacts are 
however not discussed. 



A. Setup of the NNM 

We start by replacing the underlying potential landscape, consisting of the long-ranged 
random potential due to impurities and doping atoms and the confining potential due to 
the surface of the sample, by a regular grid of nodes. This is fairly general because it is 
possible to model trajectories of arbitrary shape by states with appropriate energies at the 
saddles. On this background we replace the random potential by an effective oscillating 
potential V(x,y) = V[cos(uy) — cos(wx)] with period L := l-ajuj and amplitude V obtained 
as averages of the random potential. The potential distribution is then only topologically 
changed (e.g. deformations of contours) but edge channels and backscattering remain. 

Due to the topology of the saddle points each node connects to four chemical potentials 
and is visualized as a circle with links labeled by u\, ...,114 counter-clockwise with U\ de- 
noting the upper, right corner (see Fig. [1]). Consistency demands that only two independent 
differences of potentials exist, which represent the two components of the local electric field. 
Opposite differences are therefore equal We then define the ratio of longitudinal to 
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transversal field component, given by 

p Ex_ _ ui - u 2 

E y U\ —U4 u 2 — u 3 

In this way we construct a "transfer" equation for the chemical potentials 
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The chemical potential distribution can be calculated as a boundary value problem once the 
values for P at each node are given. External contacts supplying current are modeled by 
saddles with a pair of trajectories that point into/out of the sample (E x = 0, compare the 
Landauer-Biittiker picture and one of these two trajectories is kept fixed at the chemical 
potential of the supply contact. 

The random potential is mapped to the NNM by setting up the grid in such a way that 
in every second column the nodes are rotated clockwise by an angle of 90° and P' — 1/P is 
used for P. This reflects the topology of saddle points of the underlying oscillating effective 
potential V(x, y). As a result loops can be formed in isolated valleys or around peaks. Such 
loops have actually been observed experimentally [221 ] . 

We use the Landau levels and a self-consistent Thomas-Fermi approximation in order to 
obtain the charge density at zero temperature. Finite temperatures increase the transition 



regions as was analysed in a previous paper [23]]. The screening of the electrical bare potential 
is then obtained from the charge density simply by multiplication with a given constant 
factor of C = 50mV / (10 u cm~ 2 ) (units e = 1). Furthermore we use a constant broadening 
of 0.5meV to mimic Landau bands generated by the disorder potential. Concerning spin 
splitting a g-factor of g = 4 was used for all calculations {24]. 

We neglect effects due to feedback of the chemical potential onto the charge distribution 
and hence on the values of P. Note this amounts to the assumption that states across 
saddles describe a linear interpolation between states of two chemical potentials, that is, the 
system traverses from equilibrium to the nonequilibrium steady state adiabatically. 

In order to calculate the P-values we adopt the high-field approximation of localized 
electron states in form of contour lines of the potential, also at saddle points. In this way 
we assume a purely off-diagonal conductance tensor within the sample with a contribution 
of e 2 /h from each Landau level. Then we can immediately write P = I y /I x as a ratio of 
transversal to longitudinal component of the current. 



According to the edge channel picture we call T the probability of transmission in longi- 
tudinal direction, therefore we get I v oc R and I x oc T = 1 — R. This can be used to calculate 

n 

P from details of the Landau levels. Compare the work of Streda et. al. [25| in this respect, 
who describe the whole sample as such a single node. The current of the sample is finally 
calculated by taking the sum of all transversal potential differences in the contact nodes 
and multiply by e 2 /h. Equality of total input and output current determines the potential 
distribution uniquely. 

The present algorithm has the advantage that only the electric field has to be calculated. 
The only detail used from the Landau levels is the energy spacing of levels and the number 



of levels below Ep. Similar to other approaches [2l|, |26| we are able to calculate longitudinal 
and transversal resistance by identifying the current with the macroscopic current direction. 
The local dissipative component could be found with a formulation of tunneling and electron 
statistics in nonequilibrium together with a principle to describe the nonequilibrium steady 
state, such as, e.g., minimum entropy production, which would however make calculations 
much more demanding. Such investigations are under way and will be published elsewhere. 
In this spirit we use the expression of the tunneling transmission R mn between edge 



channels through the saddle in the presence of a magnetic field derived in Refs. 
to obtain P (Ref. [9]) 



18] and [19| 



Pr- Rmn 
= o mn _ — = exp 

1 itmn 



L 2 B 



— e 



with e := Ep — Vs the difference of the Fermi energy to the saddle energy Vs and B denoting 
the magnetic field strength. 

In the plateau region the loops in the bulk are isolated but in the transition between 
plateaus a fraction of electrons in the edge channels tunnels into the bulk, leading to a po- 
tential difference in the direction of the mainly longitudinal current and therefore dissipative 
transport. Equilibration among edge channels is allowed by simulating tunneling between 
edge channels with a chosen constant decay length for the involved states. 
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III. RESULTS 



A. Standard sample geometry 

In this paper we present a systematic study of nonideal contacts. The sample is set 
up as a typical Hall bar with two current-injecting metallic contacts and four measuring 
contacts to obtain longitudinal and transversal resistances. In the course we investigate the 
distribution of the chemical potential across the sample when varying the ratio of contacted 
cross-section of the leads as well as the length of the leads. Furthermore we placed a barrier 
across the width of the leads in vicinity of the current-injecting contacts (contact nodes), 
formed by a potential energy ridge consisting of a plateau and quadratic tails, to simulate 
non-ideal matching between the states in the contact and the 2d electron gas. The curvature 
of the saddles was chosen as a := 2hV / L 2 = 1 in all calculations of this section. 

Regarding the metallic contacts and the leads, the quantities relevant within the nonequi- 
librium steady state are the ratio r c := Wc/Wl of the width of the contact Wq to the width 
of the lead Wl and the ratio r ; := L/W, < r\ < oo, of the lead length L to its width 
W. We motivate this choice by noting that the potential distribution tends to adapt to the 
sample geometry and therefore simply scales on multiplying x and y with a common scale 
factor. 

We present an overview of transversal and longitudinal resistances against magnetic field 
in Fig. |5] by comparing samples with different r\. In Figs. |3] and [6] we show differences 
to illustrate small changes while changing r c and r^. Introducing a barrier in the current 
injecting leads is discussed in Fig. O The relative error of Hall resistances for various 
calculated setups is demonstrated in Fig. [TUJ 



1. Different metallic contacts 

We investigate two limits when applying the external potential, namely point contacts 
(r c — > 0) on the one hand and second contacting the whole cross section of the lead (r c = 1). 
Intuitively for long enough leads we expect the potential to adapt to the shape of the lead 
and the geometrical shape of the contact should be irrelevant. This has indeed been found 

n 

in experiments [27] , the rule of thumb to be sure that the measurement does not depend on 
the form of the contact is to use leads with r/ > 4. 
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In Fig. [2] we compare the potential distribution for samples with r c — > and r c = 1 
in the plateau and transition region, using r\ = 1.5 for both samples. It turns out that 
the potentials are very similar in the whole sample except near the contacts. However we 
note that in contrast to the rule of thumb the length-scale around the contact where the 
potential contours have different slopes is only a fraction of the width of the lead. In the 
plateau region this length is even smaller than in the transition region, because P m near 
the boundary forces the potential contours parallel to the longitudinal direction. 

In Fig. [3] we show the difference between resistances for point contacts and the respective 
values for r c — 1. We obtain very different magnitude of differences depending on 77 and on 
the magnetic field B. Apparently the differences are very small for 77 = 1.5 up to a field 
of about 10T. On the other hand, for 77 = the differences are much larger and one can 
clearly see AR xy tracing the structure of peaks in R xx . Furthermore it is interesting to note 
that independent of 77 and B the difference of R xx is significantly smaller than for R xy . We 
conclude that a difference between point- and line-contacts vanishes quickly with increasing 
77 and R xy is stronger influenced. 

2. Length of the leads 

We calculated five sample geometries, one with practically no leads (limit 77 — > 0) , one 
with small leads (77 = 0.5), one with medium length (77 = 1.5) and two with long leads 
(77 = 3, 77 = 5). Chemical potential distributions at selected magnetic fields are shown 
in Fig. HI To our surprise even for 77 << 4 (at least down to 77 = 1) the field contours 
adapt from the contacts to the parallel configuration in the bulk within the plateau region 
before leaving the leads, against the rule of thumb. Furthermore it seems that actually the 
lengthscale for this adjustment depends very much on the length of the leads itself. The 
adjustment is the faster the shorter the leads, again pointing to a scaling behaviour of the 
chemical potential. We can then understand this behaviour by noting that the change in 
geometry when leaving the lead, in our case the widening of the sample by the measurement 
contacts, forces the potential to adapt to the core of the sample. 

The calculated resistance differences from changing 77 restore the rule of thumb, meaning 
that resistance values converge nonlinearly with increasing 77, with 77 = 3 already close to 
the limit. In Fig. [5] we find the typical function of transversal and longitudinal resistances 
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with the magnetic field. At low field values the peaks of the longitudinal resistance clearly 
show that spin levels are not completely resolved as the 2 pairs around 2.5T and 3 — 4T 
are not yet split. We then show the difference of transversal and longitudinal resistance 
of various values of r\ to the corresponding resistances for long leads, r\ = 5 in Fig. O It 
is apparent that R xx decreases slightly with increasing r; while R xy increases at the same 
time. We also see a convergence with increasing ri, that is the differences for r\ = 1.5 are 
smaller than half the one from r/ = and only a few ohms at r\ = 3. One further notes that 
differences of R xx and R xy have opposite sign and approximately equal magnitude. Making 
the plausible assumption that the current does hardly change with ri, this points to the fact 
that the sum of transversal and longitudinal potential difference is approximately constant 
for various This can be explained qualitatively by the fact that all chemical potentials are 
evaluated by current conservation, because then the sum of potentials should stay constant 
if the current does not change. This also gives a hint, why AR XX and AR xy in Fig. [3] are 
not symmetric around zero in case of altering r c . This time we expect the current to change 
with the contact geometry. 

3. Contact barriers 

All previous calculations have been made under the assumptions that the equilibrium 
electrical potential within the leads is fairly flat. This is somehow unrealistic because the 
interface between metal and 2DEG usually develops a barrier due to charge transfer pro- 
cesses, having different potential energy. We try to mimic realistic contacts by introducing 
such a barrier in front of the metallic contacts. The height of the bare barrier was chosen 
to be 50meV at its plateau, which is reduced by screening to a value in the range 0.5meV - 
8.5meV depending on the Fermi energy. The distribution of the chemical potential is shown 
in Fig. [7] for low and high magnetic field in the plateau and transition region for the case 
of each current-injecting contact having a barrier. Furthermore we compare the scenario of 
having a barrier in only one of the leads for the transition region within the high field case. 
The most striking difference between low and high field is the larger angle of equipotential 
contours at high field relative to the longitudinal direction in the bulk for low magnetic field. 
This agrees with the increasing peak maxima of the longitudinal resistance when increasing 
the magnetic field. There are two aspects for understanding this behaviour. We note that 
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P = (transport by edge channels only) holds for all network layers except the topmost 
one. Therefore the decreasing number of Landau levels decreases the total current and con- 
sequently increases the resistances because in these layers with P = each one contributes 
the same part to the total current. On the other hand with increasing number of layers with 
P = the potential distribution is forced stronger to a plateau-like potential distribution 
even if the topmost layer has P > in the transition region. 

We compare the influence of various configurations of barriers on the longitudinal resis- 
tance in Fig. [HJ Obviously considering barriers in the measurement contacts has a larger 
effect than non-ideal current-injecting leads alone. We note that the transition starts at 
slightly lower field for barriers in all arms and the transition region is enlarged. 

In Fig. M we plotted the difference in resistances for the left lead or both leads having 
a barrier and the corresponding values for the sample without barriers. Interestingly R xx 
decreases and R xy increases when having a barrier in the leads, similar to the behaviour of 
the resistances with increasing r\. The difference between one and both leads with barrier is 
fairly small, especially for the transversal resistance. Both leads having a barrier produced 
a larger lAi?^ than one barrier, as expected. In addition we compared calculations with 
barriers of bare height lOOmeV in the current injecting leads with the analogous calculations 
for 50meV described above. The difference of resistances is significantly smaller than the 
comparison to the case of no barriers. This can be understood by the fact that the height 
after screening is not very different despite the ratio of 2 for the unscreened barriers. However 
we observe a much larger error for samples with barriers in the measurement contacts. This is 
to be expected because these barriers directly influence the potential differences measured. 
Apparently this has a larger effect than a change in total current due to barriers in the 
current-injecting leads. 

4- Dissipation 

We can learn valuable details about dissipation by inspecting the chemical potential 
distribution. We find in Fig. [2] that the equipotential contours always form loops that start 
from above a current-inducing contact and end below the same contact. Therefore, due to 
the setup that a finite current is flowing from one to the other contact it is necessary for this 
current to cross contours, meaning that the current flowing through the sample must have a 
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dissipative component. Assuming homogeneous diagonal entries to the conductivity tensor 
at saddles across the sample (which should be an adequate approximation for a constant 
equilibrium potential energy across the sample) we find that the dissipative current should be 
inversely proportional to the gradient of the potential. We conclude that entropy production 
occurs quite locally at the lower end of the left current-injecting contact and at the upper 
end of the right current-injecting contact in the geometry of our samples. This is in accord 
with measurements |28(. Especially for contacts with barriers we expect hardly any current 
crossing the barrier but rather it will follow potential contours around the barrier. This is 
always possible if the barrier is not too high, as in our simulated cases. We note that flow 
along the equipotential contours agrees with the experience that current flows along the 



edges 



291 ]. whereas the Landauer-Biittiker picture would give maximal current in the bulk 



as there are the maximum number of channels up to the Fermi energy 



30]. 



B. Nonlocal configurations and gated samples 

After investigating systematic changes at the contacted leads where current is injected 
into the sample we proceed to simulate nonlocal geometries of recent experimental investi- 
gations of non-ideal contacts with gatings. 

To set up the simulation we first generated two samples with the same geometry as used 
in the experiment. In the following the contacts are numbered by starting with 1 at the left 
boundary edge and increasing the numbers in the clockwise direction. The first sample (SI) 
had the external potential connected to contacts 1 and 2 and the voltage drop is measured 
between contacts 4 and 3. In this way we measure a longitudinal voltage difference because 
(4,3) can be connected along equipotential lines in the plateau region. The second sample 
(S2) had the external potentials connected to contacts 2 and 4 and measures the voltage 
difference between contacts 1 and 3. This therefore amounts to a transversal voltage drop. 

We ran calculations of both sample geometries by sweeping the magnetic field up to 12T. 
We used ideal samples without gating and ones where gating of some of the contacts has 
been applied in analogy to the experiment. The samples are illustrated in Fig. [11] and [12j In 
all calculations of this subsection the saddle curvature was chosen to be a := 2hV/L 2 = 0.1. 

Before discussing the potential distributions, we want to point out a general problem in 
this context. In analogy to the experiment the bulk region of the sample carries electrical 
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potential even in case it is insulating (plateau regime). This potential distribution is left 
over from the previous transition between plateaus, after which the bulk became decoupled 
from the edges due to insufficient equilibration time. Consequently the bulk potential has 
physical relevance only within the transition between plateaus. 

1. Enhanced longitudinal resistance in sample 1 

By sweeping the field we found that the sample with gating on contacts 1 and 4 shows an 
extra R xx peak centered around B = 5.7T and a significantly enlarged peak in the interval 
[8.5T, 11T], that is at the same time much broader than the respective peak of the sample 
without gates. Furthermore it forms a plateau in the middle, however with a dip around 
B = 10.33T. The plot is given in Fig. [13J It is apparent that besides the two intervals 
[5.5T, 6T] and [8.5T, 11T] the samples with and without gating show qualitatively the same 
resistance, even with slightly higher values for the ideal sample. 

To understand all these features we concentrate on the potential distribution across the 
sample, which is given in Fig. [TH for four magnetic fields within or in vicinity to the larger 
peak of the longitudinal resistance of the gated sample. At 8.00T we still have perfect edge 
channels and therefore Rl is zero. At around 8.5T dissipation sets in and Rl becomes 
nonzero. Concentrating on the distribution at 8.86T we find the reason in an equilibration 
of channels that occurs at contact 3. Around a magnetic field of 10. 3T a dip of Rl develops 
which can be attributed to equilibration among edge channels. Note the edge channel with 
the potential of contact 1 strongly equilibrates at the lower left corner with other edge 
channels which are at the potential of contact 2 and have not been equilibrated to the 
potential of contact 1 yet due to the gating. This mechanism results in less voltage drop 
between contacts 3 and 4 because the main drop occurs in the lower left corner. Finally 
around 11T the dissipation vanishes in the gated and ungated samples and we reach the next 
plateau region. Note similar equilibrations among edge channels is responsible for the non- 
ideal dissipation producing the peak starting at 8.86 T and also the peak centered around 
5.7T. 
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2. Shift of transition in transversal resistance in sample 2 

For this calculation we used sample 2 with different gating configurations. Figure [15] 
shows the transversal resistance obtained from the potential applied via contacts 1 and 3. 
There are no qualitative differences between the ungated sample and gating configurations a) 
and b). However gating all contacts leads to the interesting phenomenon of shifting the last 
plateau transition to a lower value of B by the remarkable values of about 2T. Furthermore 
for gating configuration c) a peak develops around 5.75T due to changes in edge channel 
equilibration, where all other gating configurations are within the plateau. It seems that 
also this transition starts shifted to lower magnetic fields but then falls back to the plateau 
value. This can happen when the equilibration region among edge channels changes from 
outside the range between the measurement contacts to inside. 

As an intuitive explanation for the shift of the last plateau transition it comes to mind 
that due to gating at all contacts isolated edge channels can be generated, which couple to 
only one contact. Noting that the energy of the states qualitatively follows the underlying 
potential energy surface and the Fermi energy in equilibrium is to a large extend given 
by bulk properties we arrive at the conclusion that indeed the number of occupied states is 
decreased at the gating. The plots of the distribution of chemical potentials in Fig. [161 indeed 
support this explanation. We find that at 8.3 T with and without gating the edge channels 
from contact 2 to contact 4 have the same potential, due to equilibration at contact 4 where 
the external potential is applied. On the other hand for 9.5 T the gating configuration c) 
clearly develops two edge channels between contacts 2 and 4 of different potentials, namely, 
one with the potential of contact 4 and the other of contact 2. This results in an isolated 
edge channel that cannot contribute to the current and effects the shift of the transition. 

IV. SUMMARY AND CONCLUSION 

In the present paper we summarize a series of calculations of resistances and chemical 
potential distributions for integer quantum-Hall samples with the aim to investigate the in- 
fluence of the shape and potential energy of nonideal contacts where the current enters/leaves 
the sample. In the first part of the paper we investigate the influence of the ratio of the 
contacted part of the cross-section of the leads, the length of the leads and potential barriers 
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within the leads and measurement contacts. The contours of the chemical potential more 
or less adapt to the local geometry, for instance they "flow" around barriers. From poten- 
tial contours we were able to identify the regions where dissipation will occur. Resistances 
strongly depend on the details of the sample used. We find relative errors that range up to 
a few percent in the transition region. We conclude that the Quantum Hall effect is robust 
against such changes only within the plateau regime, that is, whenever the Fermi energy lies 
between two Landau levels. 

In the second part we calculated gated quantum-Hall samples in nonlocal geometries that 
have recently been used in experiments of non-ideal contacts. Various different configurations 
of gating at the metallic contacts are considered in analogy to the experiment. We find that 
partial equilibration near a remote gated contact can lead to enhanced dissipation. On the 
other hand if all involved contacts are gated it might happen at sufficiently high magnetic 
fields that one Landau level is connected to only one of the contacts and as a consequence 
does not contribute to transport and therefore the corresponding plateau transition is shifted 
down in magnetic field. We believe that our calculations give important information for 
experimentalists, especially regarding questions of metrology. 

Recent experiments [31| used a scanned tip to apply local gating and investigated the 
change in resistance. In this way the underlying potential landscape can be locally changed. 
Using a network model that combines tunnelin g an d superconducting states the current 



distribution and hot spots have been calculated 32]. It would be an interesting extension 
of the present work to directly simulate these resistance changes, which is possible with the 
NNM. Future work is planned on this subject. 
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Figure captions 

Figure 1. Representation of tunneling at a single saddle point. The saddle is denoted by a 
circle and the lines crossing the saddle are trajectories of edge channels. The double-arrow 
denotes tunneling between the edge channels. 

Figure 2. (color online) Comparison of point contacts with contacts extending over the width 
of the current-injecting leads. Visualization of chemical potential distributions in the sample. 
The left column corresponds to the plateau region at 9T whereas for the right column shows 
the region of transition to the last plateau at a magnetic field of 11T. Positive/negative 
values of the chemical potential are shown by full/broken contours. We plotted the values 
0, ±4, ±8 with the largest /smallest value drawn in bold. The value has a separate linestyle 
for easy distinction. 

Figure 3. (color online) Difference of transversal and longitudinal resistances for point 
contacts (r c = 0) relative to values for line contacts (r c = 1), for r\ = and ri = 1.5, plotted 
against magnetic field. 

Figure 4. (color online) Visualization of chemical potential distributions across the sample 
with different current-injecting leads. The upper row shows results for the ratio length versus 
width of r\ m whereas the lower row demonstrates results for r; = 3. The magnetic field 
is set to values of the plateu (9T) and transition (11T) region. Corresponding resistances 
as functions of magnetic field are shown in Fig. [H The linestyle and values of contours is 
equivalent to Fig. [2j 

Figure 5. (color online) Longitudinal and transversal resistance versus magnetic field for 
various values r L of the length of the current injecting leads. 

Figure 6. (color online) Difference of transversal and longitudinal resistances for leads with 
various r/ relative to values for r\ = 5, plotted against magnetic field. 

Figure 7. (color online) Potential distributions for current-injecting contacts having a barrier 
in the lead. The first row shows the plateau region whereas the second demonstrates results 
for the transition region. The last row shows the case of only one barrier and the sample 
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with barriers in the leads and measurement arms. The linestyle of contours is equivalent to 
Fig. [21 however we used the values 0, ±1, ±5, ±9. 

Figure 8. (color online) Transversal resistance versus magnetic field for various configurations 
of barriers. 

Figure 9. (color online) Difference of transversal and longitudinal resistances for leads with 
and without barrier in one or both current-injecting leads, plotted against magnetic field. 
The energies denote the value of the height of the unscreened barriers. 

Figure 10. (color online) Relative Error of transversal resistances for various ri, with and 
without barrier, and the case of leads with point contacts (r c = 0), no barrier and ri = 1.5, 
compared to leads with r c = 1, no barrier and r\ = 1.5, plotted against magnetic field. 

Figure 11. (color online) This figure shows the setup of gated sample 1. The contacts are 
numbered by starting with 1 at the left boundary edge and increasing the numbers in the 
clockwise direction. The external potentials are connected via contacts 1 and 2 and the 
(longitudinal) voltage drop is measured between contacts 4 and 3. Gating potentials are 
indicated by shaded stripes along contacts. 

Figure 12. (color online) This figure shows the setup of gated sample 2 for various gating 
configurations. The contacts are numbered as in Fig. [TTJ The external potentials are con- 
nected via contacts 2 and 4 and the (Hall) voltage drop is measured between contacts 1 and 
3. Gating potentials are indicated by shaded stripes along contacts. Gating configurations 
are labelled as a) contacts 1,4, b) contacts 2,3 and c) all contacts. 

Figure 13. (color online) The longitudinal resistance Rl as a function of the applies magnetic 
field, calculated from sample 1 by dividing the voltage drop by the total current. We compare 
the ideal sample (without gatings) with the gated configuration in Fig. [TTJ 

Figure 14. (color online) Visualization of chemical potential distribution across sample 1 
as a color chart. The color label gives the respective voltage values. The magnetic field is 
indicated above each plot. 
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Figure 15. (color online) The Hall resistance Rt as a function of the applies magnetic field, 
calculated from sample 2 by dividing the voltage drop by the total current. We compare 
various possible configurations of gatings with the numbers corresponding to Fig. [TJ] and 
ideal corresponds to no gating at all. 

Figure 16. (color online) Visualization of chemical potential distribution across the sample as 
a color chart. The color label gives the respective voltage values. The upper row corresponds 
to all contacts gated while the remaining image corresponds to no gating. Magnetic fields 
are indicated above each plot. 



19 



3 % 

FIG. 1: 



B = 9.00T 




20 40 60 80 100 120 140 



B = 10.95T 




20 40 60 80 100 120 140 



B = 9.00T B = 10.95T 




20 40 60 80 100 120 140 20 40 60 80 100 120 140 



FIG. 2: 



20 



a 
pi 




R xy ,i, = 
R xy , r, = 1.5 

R„,i, = 
R.v, r, = 1.5 



10 15 20 25 

B [T] 



30 



FIG. 3: 



B = 9.00T , rj = 




10 20 30 40 50 60 70 80 



B = 10.95T , rj = 




10 20 30 40 50 60 70 80 



B = 9.00T , ij = 3 




20 40 60 80 100 120 140 160 180 



B = 10.95T , rj = 3 




20 40 60 80 100 120 140 160 180 



FIG. 4: 



21 





AR xy ,r, = 
AR xy ,i-, = 1.5 

AR xy , r, = 3 

AR M ,r, = 
AR„,i 1= 1.5 

AR.„, r, = 3 



10 15 20 25 

B [T] 



30 



FIG. 6: 



22 




B = 2.20T , both leads B = 10.95T , both leads 




B = 10.95T , left lead B = 10.95T , all arms 




20 40 60 80 100 120 140 20 40 60 80 100 120 140 



FIG. 7: 



23 



R 



xy 



no barrier 
2B, 50meV 
all arms, 50meV 



8 10 
B[T] 



12 



14 



FIG. 8: 



AR xy , IB, 50meV 
AR xy , 2B, 50meV 
AR XX , IB, 50meV 
AR XX , 2B, 50meV 
AR xy , 2B, 100 vs 50meV 
AR XX , 2B, 100 vs 50meV 
AR xy , all arms, 50meV 
,50meV 




10 
B [T] 



12 



14 



16 



FIG. 9: 



24 



6 



s 



© 




r,= 1.5, IB (vsr, = 1.5) 
r,= 1.5, 2B (vsr, = 1.5) 
r c = 0,r, = 1.5 (vsr, = 1.5) 
r, = (vs r, = 5) 
r, = 1.5 (vs r, = 5) 
r, = 3 (vs r, = 5) 



10 



15 



B[T] 



FIG. 10: 



FIG. 11: 



20 



25 



© 



© 



o 



o 



FIG. 12: 



25 



R L (Contacts 3 and 4) 



a 

O 



I/} 



14000 
12000 
10000 
8000 
6000 
4000 
2000 


-2000 




4 5 6 7 8 9 10 11 12 13 

B[T] 



FIG. 13: 



26 




27 



a 

O 



I/} 



R T (Contacts 1 and 3) 



30000 
25000 
20000 
15000 
10000 
5000 



1 1 1 1 1 1 1 r 

ideal 

1,4 gated /" "" 

2,3 gated 
all gated 



r 



I 



_i i i i i i i i_ 



4 5 6 7 8 9 10 11 12 13 

B[T] 



FIG. 15: 



28 




29 



